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1. Introduction 

Sensitivity to perturbations of quantum evolution is one of the main reasons for irreversibility 
in low dimensional quantum systems. The fidelity [1| , later dubbed Loschmidt echo (LE) [2|, 
defined as 

M(f) = \(\if \e i,H ^ h e- i,Ho/h \\i/ )\ 2 (1) 

was proposed to measure such a sensitivity. It is the overlap of an initial state | i//b) evolved 
forward up to time t with a Hamiltonian Hq, with the same state evolved backward in 
time with a perturbed Hamiltonian H^. The parameter E characterizes the strength of the 
perturbation. Equation (Q} can also be interpreted as the overlap at time t of the same state 
evolved forward in time with slightly different Hamiltonians. While the first interpretation 
gives the idea of irreversibility, the second is related to the idea of sensitivity to perturbations 
in the Hamiltonian. An important fact is that the LE, and other important related quantities 
like polarization echoes, were measured in several different experimental setups such as NMR 
(3JIU acoustic waves 0, wave guides [6], microwave billiards OH), cold atoms [9|, Bose- 
Einstein condensates [10], and quantum chaotic maps implemented with atomic spins ifTTI . 

In recent years different time and perturbation regimes of the LE have been studied in 
detail |2] [12J Q~3] [14]. We first summarize the behaviour as a function of time. For very 
short times perturbation theory gives Gaussian decay [see e.g. 031 \TE[ [TTl l. This short 
time transient is followed by an asymptotic regime. The fundamental difference between 
quantum chaotic and regular systems lies in this regime. While for the former the decay is 
exponential, the latter exhibits - on average - a power-law dependence in time [see reviews 
|[T3l[T4ll l. Finally there is a saturation given by the effective size of the Hilbert space |[T8l[T9l . 
In addition, we point out that for strongly chaotic systems and very small perturbations, if at 
Heisenberg time the decay has not reached the saturation value, then a crossover to a Gaussian 
decay is observed lfl5l[T6l . 

As a function of the perturbation strength - for chaotic systems - the decay rate of the 
LE in the exponential regime presents two different types of dependence. For small E - Fermi 
golden rule (FGR) regime - the decay rate is given by the width of the local density of states 
(LDOS), C7 LD0S ■ The LDOS is the distribution of the overlap squared connecting the set of the 
unperturbed eigenfunctions with the perturbed ones and its width is a measure of the action 
of a perturbation on the system. Then there is a crossover to a perturbation independent 
regime when <7 LD0S > A, with A the largest classical Lyapunov exponent. In this regime - 
usually called Lyapunov regime- the decay rate of the LE is given by A . This simple picture 
is supported by semiclassical results based on statistical arguments and by some numerical 
studies j2l[T3l[T4ll20ll2l1. 

Recent works j9] |22] |23] |24] EJJ |26) have nevertheless shown that the accuracy of 
the above description of the LE - which aims to provide a universal picture - is limited. 
Deviations from the perturbation independence are found usually in the form of oscillations 
around A. Moreover, there are cases where deviations are considerably large rendering the 
Lyapunov regime non-existing. 

The goal of this paper is to shed some light on the, sometimes very large, deviations from 
the Lyapunov regime. We show that, after the initial short-time transient, the behaviour of 
the LE (on average) is strongly influenced by the squared average fidelity amplitude (AFA). 
In order to do so we use the semiclassical theory known as dephasing representation (DR) 
11271 128, 29 1, and a recent result |30l , derived for local perturbations in billiards, to calculate 
the decay rate of the AFA for uniformly hyperbolic chaotic systems - in particular quantum 
maps with classical chaotic counterpart -, in the limit of fast decaying correlations - or very 
large classical Lyapunov exponents - as a function of E. 
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Our extensive numerical results show that - for quantum maps on the torus - the AFA 
imprints signatures of the perturbation to the decay of the LE not only in the FGR regime 
but also for much stronger perturbations. For a wide range of perturbation strengths, the 
semiclassical expression obtained for the decay rate of the AFA correctly reproduces the 
deviations of the LE from the - perturbation independent - Lyapunov regime. 

Aside from the dependence of the decay regimes on time and perturbation strength, 
we know that the type of initial state should be taken into account 11311 . In particular for 
semiclassical calculations the initial state must be "classically meaningful" [14] [like position 
states or narrow Gaussian states]. We show numerically that for states which are extended in 
phase space, as expected the Lyapunov regime is difficult to observe and the decay of the LE 
is dominated - at least initially - by the decay of the AFA. 

This paper is organized as follows. In Sec. [2] we introduce the average fidelity amplitude. 
In Sec. 12.11 we derive a semiclassical analytical expression the decay rate of the AFA, T^. 
Taking into account that T AFA is generally assumed to be equal to <7 LD0S , in Sec. 12.21 we 
compare both quantities. In Sec. [3]our results are tested in a paradigmatic system of quantum 
chaos, a perturbed cat map on the torus. We show numerically that for general perturbations 
(7 LD0S and T AFA can be very different (Sec. 13. U . Then, we study the decay of the LE using 
Gaussian initial states for several values of X and different perturbation types (~Sec. I3.2t . We 
show that with increasing X values the contribution of T AFA to the LE becomes more relevant. 
When the initial state is extended in phase space even for small A the initial decay rate of the 
LE is given by T AFA (Sec. I3.31 >. We summarize and expose our conclusions in Sec. [4] 



2. Average fidelity amplitude 

All the decay regimes studied previously in the literature - e.g. |2] [T2l |20l GTJ l22l . and 
reviews 117311141 - for the LE imply some kind of averaging. Typically either an average over 
perturbations or over different initial states can be done. Throughout this work we do the 
latter. We consider averaging over a number n r of uniformly distributed - random - Gaussian 
(minimum uncertainty) initial states {| Yj}Yj=v Then the average LE is given by 



M(t) = -£\Oj(t)\ 2 . 



(2) 



7=1 



where 



Oj(t) = (YrjIe^fc/V^V/) 



(3) 



is just the fidelity amplitude (FA) corresponding to state \ Yj). In this section we focus on the 
squared average fidelity amplitude [labeled AFA in Sect. [TJ 



0(t) 



7=1 



(4) 



We remark that expanding equation (@), M(t) can be expressed as n r 0(f) minus twice the 

real part of the sum of terms of type ; (f )0? (f ) with j ^ i. The average fidelity amplitude is an 
interesting quantity on its own being a measurable quantity in some of the echo experiments 
Q. It is known [e.g. in ||20l [T6] [32j [T3l [33]| to decay exponentially like 



0(f) ~exp[-r AFA f]. 



(5) 
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for chaotic systems. In general, the calculations seem to induce to the conclusion that T^ is 
proportional to the width of the LDOS ff LD0S . However in the following subsections we show 
that this result only holds either when E is small [with respect to h, see Sec. 13.11 or when the 
perturbation acts on a small [say localized] portion of the total phase space. In addition we find 
an analytical approximation T AFA which will account for the strange - strongly perturbation 
dependent - behaviour of T LE . In the particular case of quantum maps on the torus, a general 
study of the AFA and will be presented elsewhere [ 34] . 

2.1. Semiclassical calculation of the decay rate of the average fidelity amplitude 

In this section we follow the procedure introduced for billiard systems in ll30l to obtain an 
analytical expression for the decay rate of the AFA, but we extend it to a phase space setting 
(Poincare surface of section). We include a step-by-step demonstration in order for the paper 
to be as self contained as possible. One key element of this derivation is the use of the 
dephasing representation (DR) which was introduced in 11271 l28l to calculate semiclassically 
the AFA - a very thorough analysis of the DR is given in [ 29 1 . The approach is an alternative to 
that of [2] and makes the calculations more manageable. It takes advantage of the shadowing 
theorem 11291 and the initial value representation [35]. The FA in the DR is given by 

DR (t) = J dqdpW(q,p)e-' AS '^ X y h (6) 

where the action difference 

AS t (q,p,-L) = — £ f dt'V(q(t'),p(t'),t'). (7) 
Jo 

is the integral of the perturbation along the classical unperturbed orbit, and W(q,p) is the 
Wigner function of the initial state We want to average over a basis set of initial states. 
The Wigner function of the incoherent sum of any basis set is just a constant, so the average 
fidelity amplitude does not depend on the type of initial states. Therefore we have 11291 

O^Sf) = ^jdqdpe- i ^™ni h (8) 

(with Y the volume of phase space). Although the derivations in this paper can apply to any 
dynamical system in a Poincare surface of section, for simplicity the following calculation will 
be done assuming we have an abstract map on a phase space of area equal to one ("V = 1). 
Since maps are iterated in discrete steps, the time is henceforth represented by an integer n. 

Let us suppose we have a region in phase space of area a where the trajectories that pass 
through a are affected by the perturbation. The complement, with area (1 — a), is unaffected. 
Now let us divide all initial points (q,p) into sets £l n t corresponding to trajectories that have 
visited the perturbed region exactly «' times after n steps. Then equation © can be written as 

ED 

0^n)=f.OAn) (9) 

n'=0 

with 

O n i(n) = / dqdpexp[—iAS„(q,p,Y.)/h]. (10) 

Ja„, 

Due to the fact that the map is completely chaotic, the probability to hit the perturbed region 
is approximately equal to the area a. We can thus take the exponential out of the integral and 
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get ma 

O n ,(n)xi<evp[-iAS n ,(q,p,'L)/h}>a n dqdp. (11) 

Ja n , 

The integral in equation (fTTT i is just the fraction of initial conditions that visit the perturbed 
region exactly n' times after n iterations, so 



n 



dqd P = [ n , ) a"' (1 - a)"-"' . (12) 



In addition 

n' 
n=\ 

= -Z[V(q hPl ) + ... + V(q n ,,p n ,)} (13) 

where q„ and p n are just the n-th iteration of the classical map given an initial condition q, p. 
For strongly chaotic systems, q n and p n in equation (fT3l can be treated as uncorrelated-random 
variables lf30ll - at least for n small - so that the average of equation ( fTTT ) of the product of 
exponentials can be expressed as the product of averages. This assumption would be strictly 
true in the case the map was a perfect random number generator 11361 in other words, in the 
limit A — > oo. Therefore, if the average is done over a large number of initial states, comparable 
to the size of the Hilbert space, we get 

< exp[-iAS n ,(q,p,Z)/h] >n„«< evp[-iAS{q,p,Z)/h] >"', (14) 

with AS(q,p,T.) the action difference after one step. Using ( fTTT i. (1121 1. and ( TBI in equation (0 
we obtain the amplitude of the AFA 

!'A%,p,E)/fi ~\n' 



n' \ / 



= (1 -a(l- < e - iA ^.pX)/n (15) 
This exponential decay can be rewritten as 



DR 

with 



DR (n)~e- r ", (16) 



T = -ln(l - a(l- < e-mwftl* >)). (17) 
The decay rate of the squared AFA 



T DR =29*r=-21n 

AFA 



(l-a(l-<e- m( - q > p > z V n >)) 



(18) 



is then obtained by squaring the modulus of equation (IT6T ). 

The previous derivation makes equation ( fT8l 'strictly' valid in the case of local 
perturbations and uniformly hyperbolic, chaotic systems or in the case A — > °°. For local 
perturbations (05 ^ 1) a further simplification can be made by series expansion in a to obtain 
an expression equivalent to the one obtained in ll30l [37l . 

We here conjecture that equation (fT8l gives a good approximation for the decay rate of 
the AFA - for any a and A, at least for small times - based mainly on three arguments. The 
first being that in 071 a very similar derivation leading to the width of the local density of 
states is shown to hold for arbitrary values of a [see next section]. The second is that we 
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do know that for t = 1 there is no need for an a <C 1 approximation in equation (fl4l > and 
neither is there need for such approximation in the limit A — > °°. So, equation (fl4l > will be 
valid for larger times as A increases. After that there is a crossover to Lyapunov decay [which 
is observed e.g. in Fig. [2] bottom left panel]. The exact time of crossover is difficult to deter- 
mine because it depends on the form of the perturbation l34l . The last reason supporting our 
conjecture is the eloquent numerical evidence that we present further in the paper [see section 
Sec. |3| . Particularly, we show that, when the initial states are extended in phase space, e.g. 
a position state (on the torus), then even for smaller A and a 1 (global perturbation) the 
dominating decay rate - at least for initial times - is r obtained in equation ( 11 81 . This last 
point is illustrated in Sec. 13.31 



2.2. Relation between the local density of states and the fidelity amplitude 

In the bibliography of the LE |fl3l [141 it is found that the decay rate for small perturbation 
values - FGR regime- originates from the AFA and is given by the width of the LDOS (<7 LDOS )■ 
Its dependence with E is quadratic. It is the aim of this section to compare T AFA and <7 LD0S for 
any perturbation value. 

By definition the LDOS gives the distribution of the perturbed eigenstates in the 
unperturbed eigenbasis. It is well known that for small E the shape of the LDOS is a 
Lorentzian ll38l . In ll37l it is shown semiclassically that the shape - in the case of hyperbolic 
maps and billiards with a deformed boundary as perturbation - is a Lorentzian [Breit-Wigner] 
distribution 



^7) = ct-3v ; d9) 



7t((Q 2 + Y 2 

for perturbations of arbitrarily high intensity, with 

y = a(l -9? < e-^&pZ)/* >). (20) 

The result is obtained for a perturbation acting on small region of area a and then it is extended 
to any value a. The width of the LDOS can be defined as 

<^os= 2 7, (2D 

which corresponds to a distance around the mean value containing approximately 70% of the 
probability U 

We can now establish the actual relationship between the width of the LDOS and the 
AFA to the decay rate of the LE. In the case where the perturbation is local in phase space 
(a <g 1) from Eqs. (IT7l and (l20l we have 0011371 

r «a(l- < e -' A %.p)/fi ( 22) 

However, with arbitrary a the relation of equation (1231 [23ir = 2y] will not, in general, be 
true. From equation < TT~8T > it is evident that due to the logarithm, there can appear diverging 
values. 

t In fact, as seen in 1371 . the correct factor for the distance around the mean value to be exactly equal to 70% is 
1.963. Taking the factor equal to 2 gives a distance around the mean containing exactly 70.48% of the probability. 
Since this choice is in a sense arbitrary, for the sake of simplicity, we choose the latter. 
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3. Numerical results 

In this section we take a concrete system and perform numerical simulations to gauge the 
analytical results obtained for r ApA . First we check its relation with C7 LD0S . Next we do 
extensive calculations of the decay rate T LE changing various parameters like the Lyapunov 
exponent and the type of perturbation. We will also assess the influence of the type of initial 
states used for the calculations. Notably, for initial states which are extended in phase space 
(e.g. squeezed states and position states) we show that T AFA reproduces with high accuracy 
the first decay rate of the LE. We note that after this decay it is possible - depending on the 
values of both h and A - to observe a decay rate given by the Lyapunov exponent. 

For the numerical simulations we use quantum maps on the torus. Quantum maps are 
very simple systems that represent a useful tool to explore some of the main properties of 
classical and quantum chaotic systems [39 40 41 1. Furthermore, many quantum maps have 
been implemented experimentally (e.g B2ll43ll44ll45ll46l ). 

The quantization of the torus requires imposing periodic boundary conditions leading 
to a finite dimensional Hilbert space of dimension N. The associated Planck constant is 
h = \/2kN. Position basis {qi}fs l (with qi = i/n) and momentum basis {pi}fs l (with 
Pi = i/n) are related by the discrete Fourier transform (DFT). In such a setting, the evolution 
operator of a quantum map is then aN xN unitary matrix U. The operator U of the maps we 
choose to work with can be split into two operators as 

u _ e aizNT{p) e -iinNV(q) ^4) 

This form is characteristic of periodic delta kicked systems - like the standard map |47|, 
the kicked Harper map (48), or the sawtooth map. The corresponding classical map can be 
written as follows 

dV{q) \ 

dm 

1 1 dp 



(mod 1). (25) 
In particular, for this contribution we consider the perturbed cat map [see e.g. [49 1] 



p = p + aq + Kf(q) 

q = q + bp + K'g(p) (m0d1 ^' {2b) 

with a and b positive integers. For K, K 1 <C 1 this map is conservative, uniformly 
hyperbolic and completely chaotic and the largest Lyapunov exponent is A « ln((2 + a£> + 
y/ ab(4 + ab)) /2). The perturbation allows the quantum perturbed cat map to have the 
generic spectral properties of chaotic systems. The periodicities arising from the arithmetic 
symmetries of the map are eliminated for K, K' ^ ATI [49]. In particular the numerical 
results shown in this work correspond to three different perturbations: two nonlinear shears 
in momentum only 

fx (q) = 27r(cos(27r^) - cos(47r^)) (27) 

f 2 (q)=2nsm(2nq) (28) 

with g(p) — 0, and the last one composed of f\ (q) above and a nonlinear shear in position 

g(p) = 27r(sin(67rp) +cos(27T/?)). (29) 

with K' = Kin equation ( T26T ). We consider three different perturbations to illustrate more 
strongly that the analytical results obtained are robust and that r AFA correctly reproduces the 
decay rate of the LE when the system is very chaotic. 
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Figure 1. (colour online). Decay rate of the AFA compared to the width of the LDOS. In both 
panels the solid (red) line is from equation )18t and the dashed (blue) line is o"^ os of 
equation 1211 . Square symbols represent T AFA calculated from data obtained evolving n r = 10 4 
initial states for the quantum map corresponding to the classical map of equation i26\ with 
a = b = 20, and f\ (q) of equation 1271 . Hilbert space size is N = 2 18 . Circles correspond 
to the exact calculation of o" LD0S . Top panel: the perturbation is localized in the interval 
q 6 [0.25,0.35]. Bottom panel: the perturbation acts everywhere. Inset: Kf\(q); the shaded 
region shows where the perturbation was actually applied. 
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3.1. Width of the local density of states vs. decay rate of the fidelity amplitude 



The results obtained in Sec. l2.2l indicate that indeed there should be a strong correspondence 
between (T LD0S and T AFA for perturbations acting on a reduced portion of phase space of area 
a (with a <C 1). In contrast, for large alpha - in particular a 1 - the two quantities can 
differ greatly. 

To expose this property, in figure[T](top) we compare both <7 LD0S (circles), <j£JL_ (dashed 

blue line), with T^ (solid red line) for the perturbation of equation (l27l i restricted to act 
locally on a region 0.25 < q < 0.35 [see inset]. As expected the three curves are very 
close. We note that to take into account the fact that the spectrum of the cat map is periodic 
because of a compact phase space, a small correction must be introduced to the semiclassical 
approximation a sc of equation (f2Tb (see (37]). 

On the other hand, we compare af" c ' and CT T S ° with r™ (same symbols as before, and 

A LDOS LDOS AFA J 

we include T AFA - □ symbols, computed from n r = 10 4 and a = b = 20) for a perturbation 
acting on the whole phase space [figure Q] -bottom]. In the limiting case a — > 1 we have 



AFA = - 2l11 



< e -iAS(q,p)/h > 



(30) 



This leads to divergences when the argument of the logarithm approaches zero. In figure Q] 
(bottom), we see that in the global case r ApA has very high (diverging) peaks while <7 LD0S os- 
cillates with much smaller - bounded - amplitude. Only for small perturbations - E/7z < 1 - 
both quantities are approximately equal and proportional to E 2 (FGR). We have checked that 
for all the types of perturbation the results agree with the prediction. 



3.2. Decay rate of the Loschmidt echo for Gaussian initial states 

For the numerical calculations of the LE we used straightforward quantum propagation of the 
map U - equation d24l) . The propagation of this map is quite efficient due to the fact that one 
of the operators is diagonal and the other can be implemented using the fast Fourier transform. 
Given an initial state |y/b) we computed 

M(t) = \( ¥o \(ui i nu Kl r\ ¥o )\ 2 (3D 

where K\ , K 2 refer to the perturbation parameter K (= K') in the quantized version of the map 
of equation (|26| > - we take them slightly different- and E is defined by 



T.= \K 2 -Ki\. (32) 
As initial state y/Q we chose coherent - Gaussian - states 



, i y/ 4 
Woiq) = [rf&) exp 



(33) 



where | is the dispersion in q. In this section we chose ^ = 1 so the state is a symmetric 
Gaussian in phase space. The periodic boundary conditions that the torus geometry imposes 
need to be taken into account [left out of equation (l33l for the sake of simplicity]. 

We did full quantum simulations with up to dimension jV = 2 20 . We needed large 
dimension in order to have enough points, before saturation, to fit the decay rate when using 
maps with large A. In figure|2]we illustrate the different alternatives of the decay of the LE as 
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Figure 2. (colour online). LE as a function of time (n) for four different values of E for 
N = 2 20 , a = b = 4 (A ss 2.887) and using the map {26} with (q), K' = 0. The perturbation 
values are: (A) E = 4 x 1(T 7 ; (V) E = 8 x 1(T 7 ; (O) E = 1.35 X 1(T 6 ; (□) E = 1 X 1(T 5 . The 
slope of the solid red is for the corresponding perturbation - indicated with corresponding 
symbols in the inset - and the slope of the dashed blue line is A . The solid (thin) horizontal 
grey line indicates the saturation value (In AO. In the inset we show F AFA calculated from 
equation 1 1 8t and perturbation f\ of equation )27t . The points indicate the perturbation values 
taken situated at the corresponding value of T . 



a function of time in the exponential regime. We show that the decay rate can be r AFA (top- 
left, top-right) X (bottom-right), or a combination of both (bottom-left). A time evolution 
with both decay rates was also observed in [ 50 1 where, contrary to the results found here, 
the authors conclude that the larger decay rate comes from the fluctuations in the finite-time 
Lyapunov exponent. The inset shows T ApA as a function of the perturbation, and the points 
indicate the values of the perturbation for which the LE is plotted. For the following figures in 
the case where two decay rates would be visible, we take into account the one corresponding 
to smaller times [the crossover behaviour for longer times will be dealt with elsewhere l34ll l. 

Next, in figure [3] we computed the decay rate Tle, as a function of the rescaled 
perturbation T./h, for many values of a,b [see equation d26i il - i.e. different values of X and 
with perturbation f\ (q). We chose = 2 18 (we checked some of the results up to A^ 20 ) and 
averaged over 2048 initial states. For small values of X we obtain the typical result: quadratic 
dependence for small £/ ft followed by a constant equal to X (see figure [3] y), the Lyapunov 
regime. As X is increased, the Lyapunov regime becomes less visible and oscillations start 
to appear. Similar oscillations were also observed in [23, 24 25] 126). For large X, after the 
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Figure 3. (colour online). Decay rate Tle as a function of the rescaled perturbation E/S for 
the quantum perturbed cat map with N = 2 18 and for different values of a and b (different 
values of X, indicated by grey the dashed lines), (y) a = b = 1, X sa 0.96; (V) a = b = 2, 
X sa 1.76; (A) a = b = 4, X « 2.887; (□) a = b = 6, X 3.637; (O) a = b = 20, X « 5.996. 
The perturbation is /i ((<?) [equation I27H The solid red line corresponds to the r AFA , the solid 
blue line is 0^5 Averages were done over 2048 initial states. 



quadratic regime, the oscillations of T LE follow those of T AFA over a finite range. We support 
this argument with figure|4] where we show the same as in figure|3]for the two other types of 
perturbations [a combination of d2Tb with ( f29b on the top panel and (f28b on the bottom panel]. 
In both cases we see that as A increases, T LE tends to follow T AFA as a function of E. 

In both Figs.[3]and|4] there seems to be three clear regimes: if H/h < 1 then we have the 
FGR regime where T LE = <7 LD0S °= (L/h) 2 ; this regime is followed by an intermediate regime 
-which depends on A - where T LE is very influenced by r^ A - Finally for very large "L/h the 
Lyapunov regime appears. The crossover between the last two regimes needs further study, 
which would include a general understanding of the AFA as a function of the perturbation 
strength for moderate A values and also understanding the contribution of the "cross" terms 



in the expansion of 



0(t) 



3.3. Extended initial states 

In the previous section the numerical results show clearly that the exponential decay of the 
LE for large A depends crucially on F AEA . In this section we show this can also be the 
case for small A. This happens when the type of initial state chosen has a large linear 
extension. Already at the beginnings of LE research, the universality of the perturbation 
independent regimes was questioned and shown to be deeply dependent on the type of initial 
state [3J]. Here we show for classically meaningful initial states [ 14], like position states, and 
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Figure 4. (colour online). Decay rate Tle as a function of the rescaled perturbation T./H for the 
quantum perturbed cat map with N = 2 18 , and for two different perturbations [see Eqs. {26), 
(27), (28),{29) ]: (top) Kfi (q), K'g(p) and (bottom) Kf 2 {q). The different symbols represent 
different values of a and b (different values of A, indicated by the dashed lines), (top): (A) 
a = b = 6,X^ 3.637; (O) a = b = 20, A 5.996; (□) a = b = 50, X a 7.824. (bottom): (v) 
a = b = 1, X « 0.96; (O) a = b = 2, X ss 1.76; (A) a = 6 = 4, A « 2.887; (□) a = 6 = 6, 
A « 3.637; (O) a = b = 20, A « 5.996. The solid red line corresponds to T AFA , the dashed 
blue line is <T S ^ Averages were done over 2048 initial states (Gaussian). 




Figure 5. (colour online). The white circles indicate the decay rate Tle as a function of the 
rescaled perturbation E/fi for the quantum perturbed cat map with a = b= 1, N = 2 12 . Average 
is done over ~ 10 3 initial position states. The red line corresponds to r^J! , while the dashed 
blue line indicates o~ s ^ . The dashed horizontal (black) line indicates A ks 0.96. 



for squeezed Gaussian states that the decay rate of the LE is given by r even for small 
values of A . 

We first choose position states. On the quantized torus position states can be pictured as 
strips of length one and width h/(2n). In figure [5] we computed T LE where the initial states 
used were random position states. The map is the perturbed cat map described at the beginning 
of Sec. [3] with a = b = 1 corresponding to A « 0.96. We observe that T LE obtained coincides 

almost perfectly with r ARA . For the sake of clarity and contrast we also included the plot for 
°ldos ( s °lid blue line). We remark that in the case of extended initial states our conjecture 
that T AFA of equation (IT8b gives the correct decay rate of the LE is fulfilled for initial times, 
for arbitrary values of A and a. We would like to also point out that after a transient time, 
if the Hilbert space is large enough, then a second regime with decay rate given by A might 
also be observed. However for large A the Lyapunov regime becomes increasingly difficult to 
observe because of the saturation at 1 /N. 

Another way to illustrate the dependence on the initial state is by varying E, in 
equation d33l . We then get squeezed coherent states. If % = 1 we have the case of Sect. 13.21 
For 4 < 1 the state squeezes vertically and stretches horizontally, approaching a momentum 
state for % <C 1 ; for S, > 1 the state squeezes horizontally and stretches vertically, approaching 
a position state for % 3> 1 . In figure|6]we show Tle as a function of the squeezing parameter | 
-starting from a momentum state [£, <C 1], passing through a circular coherent state [£, — 1], 
over to a position state [| > 1]. We plotted this for a = b = 4 (A &s 2.887), N = 2 14 and 
10 3 initial states in the case of perturbation f\(q) [equation ( 1271 )1 . We took two different 
values of the rescaled perturbation: E/7z = 2.467 (□) corresponds in figure [3] to a region 
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Figure 6. Tle as a function of the squeezing parameter ^ for N = 2 14 , a = b = 4 and two 
different values of £/S: (□) E/fi = 2.467; (O) Y./h = 46.55. The average was done over 
n r = 10 3 inital states. The dashed line indicates X = 2.887. The solid horizontal (red) lines 
indicate the corresponding value of r^J! for each perturbation. The panels on the top show the 
Husimi function of the initial states corresponding selected values of the squeezing parameter 
(a) 0.01; (b) 0.1; (c) 1; (d) 10; (e) 100. The axes range of both (q,p) in the Husimi function 
plots is [0, 1]. 



where T LE follows T afa and thus the value of the decay rate remains approximately constant; 
T./H = 46.55 (O) corresponds to a region where Tle has already settled around A. In this case 
we observe a transition from a decay rate dominated by the AFA for £, < 0.1, £, > 10, i.e. 
where the coherent state has been stretched in either in q or in p by a factor of the order 10. 

4. Conclusions 

The Loschmidt echo is a quantity which was introduced as a measure of the instability and 
irreversibility of a system under an external perturbation. The Lyapunov regime of the LE 
describes how for some parameter values the decay becomes independent of the perturbation 
and the decay rate is the Lyapunov exponent of the classical system. It is a well known 
theoretical prediction that can be seen in simulations [e.g. for the kicked top 0121 . in the 
stadium billiard BP . or a two-dimensional Lorentz gas lED l but which for some systems 
that exhibit generic chaotic behaviour - like quantum maps - can be truly difficult to find due 
to very strong deviations J5] |22] [23] |24] [25] [26) . In this work we have thus set out to try to 
understand this elusive nature of the Lyapunov regime in quantum maps. For these systems 
we show that the AFA plays a significant role and we dare conjecture is a greater defining 
feature of a system. A similar conclusions could be inferred from previous works on billiards 
with perturbed walls Il30ll3~7l . We think the extension of our results to more realistic models, 
which fulfill the uniformly hyperbolic and chaotic condition, should be straightforward. 

We have derived a semiclassical expression r AFA for the decay rate of the AFA in the case 
of highly chaotic maps (large X limit). Using extensive numerical simulations in quantum 
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maps we have shown that, the decay rate of the LE is very well reproduced by T afa for some 
range of perturbation values. As a consequence, there are three regimes depending on the 
perturbation. For small L/h there is the well known FGR regime. Then there is an oscillating 
behaviour which can be explained with the decay rate of the AFA and finally for very large 
JL/h there is the Lyapunov regime. We point out that the decay rate obtained for the AFA 
correctly explains the oscillations that are observed numerically in this paper but which had 
been also observed previously 11231 l24l l25l l26l . The complete understanding of the crossover 
between the last two regimes is left for further studies [34|. 

We have also compared r ApA to the width of the LDOS and indicated their fundamental 
differences when the perturbation is not localized. Moreover, we have shown that the influence 
of r ApA in the decay of the LE is enhanced when the initial states are extended in phase space. 

We remark that the reasons we find here to the absence of a Lyapunov regime are 
fundamentally different to those found in ll22l where the fluctuations due to different ways 
of averaging are shown to lead to a double-exponential decay of the LE. 
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